Experimental evidence of shock mitigation in a Hertzian tapered chain 
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We present an experimental study of the mechanical impulse propagation through a horizontal 
alignment of elastic spheres of progressively decreasing diameter cj) n , namely a tapered chain. Ex- 
perimentally, the diameters of spheres which interact via the Hertz potential are selected to keep 
as close as possible to an exponential decrease, cj> n +i = (1 — q)<f>n, where the experimental tapering 
factor is either qi ~ 5.60 % or ~ 8.27 %. In agreement with recent numerical results, an impulse 
initiated in a monodisperse chain (a chain of identical beads) propagates without shape changes, 
and progressively transfer its energy and momentum to a propagating tail when it further travels 
in a tapered chain. As a result, the front pulse of this wave decreases in amplitude and accelerates. 
Both effects are satisfactorily described by the hard spheres approximation, and basically, the shock 
mitigation is due to partial transmissions, from one bead to the next, of momentum and energy of 
the front pulse. In addition when small dissipation is included, a better agreement with experiments 
is found. A close analysis of the loading part of the experimental pulses demonstrates that the front 
wave adopts itself a self similar solution as it propagates in the tapered chain. Finally, our results 
corroborate the capability of these chains to thermalize propagating impulses and thereby act as 
shock absorbing devices. 

PACS numbers: 83.80.Fg, 05.45.-a, 62.50.+p 



I. INTRODUCTION 

Granular materials are useful for shock protection, and 
cast iron shots are a well known example of system effec- 
tively used in damping of contact explosive loading 0. 
Recent numerical results have demonstrated that new 
features can still be developed, such as impulse thermal- 
ization 0, 0, 01- or shock confinement [|| in one dimen- 
sional granular media. Indeed, the propagation of an im- 
pulse through a chain of monodisperse elastic spheres ex- 
hibits interesting nonlinear physics 1] . When the spheres 
in the chain barely touch one another, the energy of an 
impulse initiated at one end of the chain propagates as 
a solitary wave of well defined velocity and width, as de- 
scribed by pioneering work of Nesterenko Ju El pl E j , E3 , 
and as confirmed later by other authors [ll], [13, [13, uM ■ 
When small dissipative effects such as viscoelastic resti- 
tution or solid friction (mainly thwarted rotation between 
beads and friction between beads and walls of the setup) 
are taken into account the solitary wave only attenu- 
ates in amplitude and spreads 0, [^. Under static 
loadin g, t he chain still supports solitary wave solutions 
0, ITU Il3l| , in addition to dispersive linear sound waves 
of weak amplitudes. In the presence of a gradient of 
static loading (e.g. a vertical chain under gravity), it 
has been proved that propagation of an impulse is con- 
trolled by dispersion, and extra dispersion introduces a 
coupling between quasisolitary and oscillatory propagat- 
ing modes [lTJ H, EH • However, dispersion does not 
allow the possibility of distributing the energy of a soli- 
tary wave uniformly, neither throughout a finite length 
vertical chain nor m a horizontal linear chain [l3| . 
In regard to energy propagation in the absence of disper- 



sion, Poschel and Brillantov [ljj considered a one dimen- 
sional chain of spheres with a restitution coefficient that 
was velocity independent. By keeping constant the size 
of the beads at both end of the chain they showed that 
the energy transmission is optimal when the mass distri- 
bution is an exponentially decreasing function. Similarly, 
Nesterenko et al have shown 0,0,13 that the momentum 
of a pulse is totally transmitted when it passes through 
a sharp decrease of bead's size, and is partially reflected 
in the reverse case. 

Here, we investigate experimentally the effect of shock 
mitigation occurring in an unloaded chain of exponen- 
tially decreasing bead size. This effect, recently predicted 
by Sen et al U0, 0, , is due to the progressive delay 
on the transmission of a fraction of the energy of the 
pulse when propagating in a sequence of increasing stiff- 
ness contacts defined by the sequence of smaller beads. 
Consistently, excepting the reflection of the pulse at the 
rigid end of the chain, there is no detectable reflected 
wave propagating backward while the pulse propagates 
down toward the smaller extremity of the tapered chain. 
Experimental evidence of this effect have been observed 
by Nakagawa et al pfl | , when an impulse is introduced in 
a tapered chain in which the beads are initially barely in 
contact and the chain is let to expand after the impulse 
propagation, the last beads being initially at some dis- 
tance of a wall sensor. This configuration is suitable to 
measure the speed of beads by imaging techniques and 
the force due to the ejection of the last bead colliding 
with a wall sensor. In our experiment the beads always 
keep in contact which allows us to adopt a complemen- 
tary point of view and focus on the shape and speed of 
the wave instead of that of beads. The main advantage 
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is that the wave picture is better realized and energy 
transfer from the front pulse to a propagating tail is fully 
visualized. This article is organized as follow. In sec- 
tion II, we describe our main experimental considerations 
and present force measurements acquired in two differ- 
ent chains, whose tapering factor is either q\ ~ 5.60 % 
or q2 ~ 8.27 %. Similarly, the acceleration of the front 
pulse is characterized as it propagates toward the de- 
creasing end of the chain. Effects of misalignments and 
tapering factor variations are investigated as well. The 
section III is devoted to brief summary of the main the- 
oretical aspects and discussion. We show that the ballis- 
tic approximation in conjunction with the effect of small 
restitution coefficient is well suited to capture the accel- 
eration and the decrease of force amplitude of a pulse 
traveling through the tapered chain. 



II. EXPERIMENTAL OBSERVATIONS 

Beads are high carbon chrome hardened A I SI 52100 
steel roll bearing. Density is p — 7780 kg/m 3 , Young's 
modulus is Y — 203 ± 4 GPa, Poisson ratio is v = 0.3, 
and yield stress is ay — 2 GPa, We have carefully 
checked that bead's deformation always keeps elastic be- 
havior in all our experiments. Assuming that the contact 
surface is a disk of area A = tt(ORF) 2 ^ 3 , we estimate 
that yield stress occurs for a compression force greater or 
equal to Fy ~ 470 N, which corresponds to an overlap 
by ~ 11 /im, for beads of 26mm diameter. This threshold 
is greater for smaller beads. 

A solitary wave is initiated in a monodisperse chain 
made of uq — 16 beads whose diameter is <p$ = 26 mm, 
by applying a short impact with constant velocity (over 
all experiments) of a smaller bead {(pi — 8 mm) at one 
of its extremity ^4|. This impulse then reaches a ta- 
pered chain composed of beads whose diameters 4> n pro- 
gressively decrease. More precisely, two tapered chains 
are implemented, whose intermediate diameters, <p n , are 
given in the table [I] Tapering factor is determined as 
q = 1 — ((0n/0o)"~™ o )n>ri o - As everywhere in this doc- 
ument, errorbars refer to the unbiased sample standard 
deviation (confidence interval, CI ~ 68 %). The first 
chain provides a tapering factor of q\ = (5.60 ± 0.67) %, 
and the second one provides a tapering factor of qi = 
(8.27 ± 0.31) %. 
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TABLE I: Bead's diameters. 

Beads of the monodisperse chain are aligned on a rigid 
Plexiglas track over which beads are allowed to roll with 
little friction. Beads of the tapered chain are located 



on a special stepped track composed of a collection of 
short tracks of length nearly the beads diameter, designed 
in such a way that the tapered chain is automatically 
aligned on the axis of the monodisperse chain. The whole 
rigid structure serving for beads alignment is fixed on 
an optical table, which guarantees a misalignment, i.e., 
off center error, lower than 2 % between the smallest 
consecutive beads, and lower than 0.1 % for the biggest 
beads considered here. It appears that repeatability of 
measurements is not so much affected by misalignments if 
the beads are not allowed to move laterally. This effect is 
quantified further in the text by introducing a controlled 
off center displacement of the beads. 




Frequency [Hz 



FIG. 1: (a)-(g) Force as a function of time (arbitrary origin). 
Evolution of a pulse initiated in a monodisperse chain of no = 
16 beads, shown in (a), as it propagates in the tapered chain of 
ni beads with tapering factor q — q\. (b)-(g) Actual forces felt 
at the right contact of beads #18, 20, 22, 24, 26, and 28 (see 
table HJ. (h) Power spectrum, with arbitrary unit, indicating 
a comparison between incoming and outcoming pulses in the 
tapered chain. 

A piezoelectric dynamic impulse sensor (PCB 208 All 
with sensitivity 112.40 mV/N) located at the end of the 
chain provides the force at the rigid end. This sensor 
has a flat cap made of the same material as the beads. 
Forces inside the monodisperse chain are monitored by 
a flat dynamic impulse sensor (PCB 200B02 with sen- 
sitivity 11.24 mV/N) that is inserted inside one of the 
beads, cut in two parts. The total mass of the bead sen- 
sor system has been compensated to match the mass of 
an original bead. This system allows achieving non intru- 
sive force measurement, by preserving both contact and 
inertial properties of the bead-sensor system. The stiff- 
ness of the sensor k s = 1.9 kN//im being greater than 
the stiffness of the Hertzian contact (k s ^> kn oc kS 1 ^ 2 , 
see section III for definitions of n and S) , the coupling be- 
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FIG. 2: Same as Fig. Q for a pulse initiated in a monodis- 
perse chain of no = 16 beads, shown in (a), and propagating 
through the tapered chain of n.2 beads with tapering factor 
q — q2- (b)-(g) Actual forces felt at the right contact of beads 
#18, 20, 22, 24, 26, and 28 (see table Q. (h) Power spec- 
trum, with arbitrary unit, indicating a comparison between 
incoming and outcoming pulses in the tapered chain. 



tween the chain and the sensor is consequently negligible. 
Additional details on how to relate the force F s detected 
by the sensor with the actual force at the contact be- 
tween two beads are given in Ref. fyj. Signals from 
sensors are amplified by a conditioner (PCB 482A16), 
recorded by a two channels numeric oscilloscope ( Tek- 
tronix TDS340), and transferred to a computer. Mea- 
surements have been repeated nine times and averaged 
to minimize errors, and acquisitions are achieved as fol- 
lows. In the monodisperse chain, the sensor at the end 
of the chain serves as a trigger, and the force is recorded 
at various positions by the sensor in the chain. In the ta- 
pered chain, force is recorded by the sensor at the end of 
the chain, while the sensor in the monodisperse chain is 
let at a given position to trig the acquisition. In addition, 
we know from numerical simulations 01 that under per- 
fectly elastic collisions, forces recorded in the chain and 
at the end of the chain are related in amplitude and du- 
ration: F™ atl ~ 1.94 x F^ am and t wM ~ 1.09 x Tchain , 
with relative error less than 5 %. 

In Fig. n the evolution of an unaveraged pulse, initi- 
ated in a no — 16 beads long monodisperse chain, as 
it propagates in the first tapered chain (q = qi) made 
with n\ — 14 beads of decreasing diameters as described 
above. Practically, the evolution of the pulse is measured 
at the end of a tapered chain containing an increasing 
number of beads (0 < ni < 14). It appears that a tail is 
formed behind a front pulse and the wave, initially sym- 
metric, becomes more and more asymmetric. In Fig. |2 



is the evolution of a pulse, initiated in a uq — 16 beads 
monodisperse chain, as it propagates in the second ta- 
pered chain (q = q 2 ) made with n 2 — 12 beads of decreas- 
ing diameters as described above. The tail formed behind 
a front pulse of the wave is in this case even more pro- 
nounced. This effect is exactly the one predicted by Sen 
et al For tapered chains (g ^ 0), the pulse becomes 
distributed throughout the chain as it propagates. Since 
impulsion is transmitted down the chain at the contact 
surface, impedance mismatches are introduced by taper- 
ing the diameters, every mismatch limiting transmission 
between the grains. By comparing Fig. ^to Fig. [21 it 
can be observed that the effect of pulse spreading and 
maximum force decrease become more dramatic as the 
tapering factor is increased. Fig.QJi and Fig.|2fi illustrate 
typical power spectrum of the incident solitary wave, as 
well as the corresponding waves reaching the end of the 
chain for the values of tapering factor used here. A clear 
cascade of energy is observed, from low to high frequen- 
cies, as the pulse propagates through the tapered chain. 
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FIG. 3: Force as a function of time, measured at the end of 
the entire chain of of tapered factor q = 52- Vertical displace- 
ment e of the bead under scope is indicated on the figure, 
(a) Effect of a controlled vertical misalignment introduced at 
the bead number 022 = 16 mm close to the middle of the ta- 
pered chain, (b) Effect of similar vertical misalignment which 
is this time located at the bead in contact with the wall sen- 
sor, </>28 = 9 mm. (c) The force as a function of time, mea- 
sured at the end of a partial tapered chain (q — qi) made of 
beads 01 to 24. dq/q — +1 corresponds to the same signal 
when the bead number 18 is replaced by a bead of diame- 
ter 0i8 = 0i9 = 23.0 mm. Thus, two consecutive beads have 
identical diameter. Similarly, for dq/q = —1, the bead number 
19 is replaced by a bead of diameter 19 = 0is = 23.0 mm. 

In order to check the robustness of our experimen- 
tal configuration, we have explored how sensitive is the 
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shock mitigation effect to small misalignment of beads. 
To quantify the effect of beads misalignments we have 
performed two test experiments. In a first experiment, 
we check the chain sensitivity to misalignment occurring 
in the middle of the tapered chain. Second, we impose 
a well controlled vertical displacement of the bead in 
contact with the wall sensor and observe how the pulse 
is modified. As seen in Fig. |3Ji, as the misalignment 
in the middle of the tapered chain increases, the sig- 
nal decreases slightly in amplitude and the oscillations 
of the tail progressively decrease. When looking care- 
fully, one can see that the more pronounced peak of the 
tail (t — T ~ 0.38 ms) is progressively decreasing in am- 
plitude, suggesting the existence of a correlation between 
the oscillating tail and a given contact at the chain. Sim- 
ilarly, as observed from Fig. a misalignment up to 
10 % at the end chain does not introduce dramatic effect 
on the wave. The amplitude is slightly modified and only 
some effects take place at the tail of the wave. However, 
in both cases, the front pulse velocity is nearly unchanged 
at the experimental resolution achieved here. 

Although the tapering factor was not kept exactly con- 
stant through the chain in the above description, the ef- 
fect of shock mitigation predicted by Sen et al Q is well 
captured in our experimental results. To better char- 
acterize the effect of local variation of tapered factor, 
we have tested a short tapered chain in which we re- 
placed one of the bead by a smaller one, in such a way 
that two consecutive beads have equal diameter. Thus, 
a local variation of q, dq/q = +1 is produced. Simi- 
larly, to investigate the effect of a negative variation of q, 
dq/q = — 1, we replace a bead by one equal to the bigger 
neighbor. As seen in Fig. two effects are visualized. 
First, for dq/q = +1, the front pulse is slightly steeper 
and narrower. In contrast, for dq/q — —1, a smoother 
and slightly wider front pulse is produced. These two 
effects are consistent with results presented in Fig. ^ and 

Fig.ua 



III. DISCUSSION 

The physical behavior of solitary waves in chains of 
N identical beads is summarized below. Under elastic 
deformation, the energy stored at the contact between 
two elastic bodies submitted to an axial compression cor- 
responds to the Hertz potential j^, Uh = (2/5)k(5 5 / 2 , 
where 8 is the overlap deformation between bodies, k _1 = 
(e+e^iR-^+R'- 1 ) 1 / 2 , 6 = 3(l-i/ 2 )/(4Y~), and R and R' 
are radii of curvature at the contact. Y and v are Young's 
Modulus and Poisson's ratio respectively. Since the force 
felt at the interface is the derivative of the potential with 
respect to S, (Fh — dslln = nS 3 ^ 2 ), the dynamics of the 
chain of beads is described by the following system of N 
coupled nonlinear equations, 



where overdots denote time derivatives, m is the mass, 
u n is the displacement of the center of mass of bead n, 
and the label + on the brackets indicates that the Hertz 
force is zero when the beads loose contact. Under the 
long-wavelength approximation, X R (where A is the 
characteristic wavelength of the perturbation), the con- 
tinuum limit of Eq. ^ can be obtained by replacing the 
discrete function u n ±\{t) by the Taylor expansion of the 
continuous function u(x ± 2R,t). Keeping terms of up 
to the fourth order spatial derivatives, Eq. pleads to the 
equation for the strain tp = —d x u > 0, 

4, ~ c 2 d xx [^ 2 + (2/5)i?V /4 d ra (V 5/4 )], (2) 

which admits an exact periodic solution in the form of 
a travelling wave, jM£ = X — vt), with speed v. This 
solution is m H I! CO JSl : 



(3) 



= (5/4)>/c) 4 cos 4 [£/(i?V10)] 



where c = (2R) 5 '' i (n/m) 1 ' 2 . Although this solution only 
satisfies the truncated Eq. ^ it is well established that 
one hump (— 7f/2 < £/(i?vT0) < n/2) of this periodic 
function represents a solitary wave solution 0, 0, uM ■ 
Finally, approximating spatial derivative by finite differ- 
ence, the strain in the chain reads tp ~ S/(2R). The 
force felt at the interface, F ~ k{2R^) 3 / 2 , and the wave 
velocity v become 



F m cos 6 



(5^) (pf) 5 " \5wp) 



(4) 
(5) 



Ki-i 



3/2 
+ 



3/2 



(1) 



The second expression for the wave velocity v is ob- 
tained by relating the maximum force F m to the maxi- 
mum bead's velocity V m [J, 0, HI 123 ■ Experimentally, the 
maximum force F m , the time of flight T, and the duration 
t, are obtained here from fitting the loading part (— r < 
t — T < 0) of the measured force as a function of time to 
the Nesterenko solution F(t) = F m cos 6 [(f - T)/r}. 

It is worth mentioning that Rosenau et al, while study- 
ing the role of nonlinear dispersion in pattern formation, 
derived a Korteweg-de Vries like equation |2^, |27| . 
which is similar to the continuum limit of Eq.^described 
above. Among other solutions such as kinks, solitons 
with infinite slope, periodic waves, dark solitons with 
cusp, and other non analytical solutions |27j . they found 
a solution of the form of Eq. [3] with a compact support 
(thus named compactons): their solution is strictly zero 
outside a finite range, i.e. for |£/(.R\/T0)| > n/2 p7| . 

In the following, we discuss what features of the front 
pulse, travelling in the tapered chain, can be captured 
using the simple monodisperse solution of solitary waves 
in conjunction with a ballistic approximation. Indeed, 
Doney et al 4] proposed recently a ballistic approxima- 
tion, with an heuristic term rendering energy dissipation, 
to model elastic collision of individual beads. Taking into 
account a classical restitution coefficient, e, the ratio of 
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consecutive bead's velocity is, 

yn+l 



2c 



yn 



1 + (1 - qf 



(6) 



The restitution coefficient is obtained from force mea- 
surements acquired in the monodisperse chain (q = 
0), by using that e = (F« +1 /F£) 5 / 6 = V£ +1 /V£. We 
thus estimate that e ~ 0.985. The theoretical maximum 
of the force in the tapered chain is then obtained from 
Eq.0 



F>' 



2c 



16/5 



(7) 



and the predicted speed of the quasisolitary wave (the 
front pulse), in the ballistic approximation, reads, 



'Jn+l 



2c 



T + (1-9)3 



11/5 



(8) 



In writing Eq. [7| and Eq. [H] it is assumed that the 
front pulse adapts itself to the self similar solution at 
every impedance mismatch. Therefore, this approxima- 
tion is suitable for small tapering factor, since for slow 
variation of bead size, quasisolitary wave adjustments 
are expected to be small and to take place at distance 
much longer than its characteristic width. This can be 
shown also by comparing the time scale of the wave, 
T W ave ~ 2R/v, to the time scale of the loading during 
a collision, t co u ~ S m /V m , (use Eq. to estimate the 
ratio v/V m ). These characteristic times are equal in a 
monodisperse chain and almost equal for small tapering 
factor, demonstrating that the loading of the contact is 
the process governing the main features of quasisolitary 
wave formation. In Fig.^J when the force is normalized to 
the maximum force, F m , and the time scale is selected to 
be the measured duration r, all data presented in Fig. ^ 
and Fig. [21 show a collapse in the rising part of the front 
pulse. Thus, this result provides strong evidence that the 
main front of the pulse follows a self similar solution. 




Normalized Time (t-T) A 



FIG. 4: The collapse of all data presented in Fig.Qand Fig. [2] 
when the force is normalized to the measured maximum force, 
F m , and the time scale is (t — T)/t, where T and r are the 
measured time of flight and pulse duration, respectively. 

Based on the results above, we discuss in further de- 
tail the main features of the quasisolitary wave, namely, 



its maximum force, its local speed and the momentum 
transfer from the main peak to its tail. As a reference, 
for both tapering factors, Fig. [SJil and Fig. |SJi2 indicate 
the size of beads as a function of position. In Fig. [3)1 and 
Fig. [5^2 the maximum force felt by the rigid wall. An im- 
portant decrease in the force is observed in both cases as 
the pulse propagates in the tapered chain. Interestingly, 
such a decrease can not be attributed to the dissipation 
only. Indeed, dashed line in Fig.[SJ»l and Fig. 0>2 indicate 
that predictions from Eq.0 with q = 0, overestimate the 
force. However, when tapering factor is introduced, the 
experimental data follow reasonably well the predictions 
from Eq. [7| 
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FIG. 5: . Characteristics of the chain and of the pulse trav- 
eling in both tapered chains (the left column is for q = qi 
and the right one is for q = 52) as a function of the position 
in the chain. (al)-(a2): the beads diameter used in the ex- 
periments compared with the fit presented in the text, (bl)- 
(b2): maximum compression measured in the chain, (cl)- 
(c2): duration of the pulse contrasted to predictions from dis- 
sipative ballistic approximation including or not the tapering 
effect. (dl)-(d2): front and tail linear momentum, respec- 
tively Q F = J^ T F(t)dt and Q T = ~ F(t)dt. Here the 
solid lines are guides for the eye. 



In Fig. [SjH and Fig. [3i2 we represent the duration 
t of the pulses as the bead diameter decreases. From 
the theory of monodisperse chain schematized above, it 
is known that the quasisolitary wave duration r, ac- 
cording to Eq. is T n = (^„/2)VT0/t) n . It has been 
checked that r n , measured at every bead, agrees with 
the monodisperse-ballistic approximation in a reason- 
able manner, for small tapering factor, see solid lines on 
Fig-IStl and Fig. [SJi2. This result is not obvious since to 
calculate the solid line on Fig.[SJ;l and Fig. 0:2, we have 
used the calculated value of the quasisolitary wave speed 
at every bead from Eq. [SJ 
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On the other hand, the impulse of the front is Qf = 
J^ +T F(t)dt and the impulse of the trail is Qt — 

$t+t F{t)dt- The principal effect of the tapered chain 
is then shown in Fig. 03^1 and Fig. 03?2, in which a no- 
ticeable part of the impulse is transferred from the front 
pulse to the tail. From comparison to Fig.[SJ> and Fig.jSJ;, 
it appears that the diminution of the front impulse is 
principally due to both dramatic decreases of the maxi- 
mum force and of the duration of the front pulse, as it 
propagates down the tapered chain. Note that the total 
momentum is conserved in the tapered chain. This has 
been checked more accurately using lower frequency res- 
olution, which allowed us to fully record the pulse tail 
and thus verify that (Qf + Qt = cte). In fact, the sig- 
nal being acquired on a fixed set of points, for the rest 
of our measurements, we have preferred to keep a high 
sampling frequency, which provides high resolution for 
the front pulse but avoids to record the very end of its 
tail. 
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q = 0. The decrease observed for both tapered factors 
can not be attributed to the restitution coefficient effect 
only. The effect of tapering factor is better observed in 
Fig. |SJd and Fig. [f)J>2 which provide a close comparison 
of time difference, T — T(e;q = 0), where T(e:q = 0) is 
the prediction that only includes the restitution coeffi- 
cient effect. It is seen that the measured time of flight is 
shorter than the prediction from ballistic approximation 
which slightly overestimates the time of flight between 
two consecutive beads. 

For the sake of comparison, we finally contrast theo- 
retical velocity from Eq. |S] to local velocity v% — (S n — 
i)/(T n — T n -i). Experimentally, the local velocity, es- 
timated from time of flight differences, naturally provides 
very large relative errors. We thus choose to present an 
estimate obtained by smoothing the local time of flight 
fluctuations. This is done as follow. In accordance to 
Eq. [SI we map the local velocity to the behavior 
v n = Q n ~ n Vn 1 t° find equivalent coefficients of propor- 
tionality Q for both tapered chains. Results of the fits 
are presented in Fig. [Ul and Fig. 0:2. The inflection 
of the velocity in the tapered part of the chain, render- 
ing the acceleration of the pulse, is clearly demonstrated 
in both cases. Nevertheless, experimental fits slightly 
overestimate theoretical predictions, and a closer com- 
parison is necessary. Experimentally, we find that the 
coefficient of proportionality is Q\ = 1.025±0.016 for the 
first chain (q = gi), and Q 2 = 1.049 ± 0.044 for the sec- 
ond chain (q = (72)- These coefficients must be compared 
to their theoretical predictions, respectively Q\ h — 1.014 
and Q l 2 — 1.022. We thus check that both theoretical 
values lie inside their respective experimental errorbars. 
Moreover these experimental coefficients only overesti- 
mate their theoretical counterparts by 1.1 % for the first 
chain, and by 2.7 % for the second one. This last satis- 
factory agreement, between theoretical and experimental 
coefficients of proportionality Q in both chains, is not sur- 
prising and must be balanced, since Q is weakly sensitive 
to tapering factor variation 5q (see Eq. |S]with q << 1): 
SQ/Q = (3/10)<5g. 



FIG. 6: Time of flight and velocity measurements, in both 
tapered chains, as a function of the position in the chain. 
(al)-(a2): flight time T. (bl)-(b2): difference between the 
measured flight time and the theoretical flight of time ob- 
tained when only the restitution coefficient is taken into ac- 
count (e < 1 and q = 0). (cl)-(c2): pulse velocity. 

We compare below the measured to the predicted 
time of flight of the front pulse when reaching the sens- 
ing wall. The relation T n = To + J2i (<t>m/ v m) = 
To + (4>o/v\) \Z- € -i/5 helps to find v\ and Tq, by fit- 
ting to measurements done in a chain of equal beads 
(n < no). The theoretical velocity is then found by us- 
ing Eq. |SJ and finally, a prediction of the time of flight 
in the tapered chain is obtained by using the relation 
T n = T n _i + (4> n /v n ). In Fig. EJil and Fig. EJi2 compar- 
isons with experimental data are presented in the case of 



IV. CONCLUSION 

In conclusion, we have shown that a tapered chain is 
useful for shock attenuation and energy transfer from 
low to high frequencies, as previously stated by Sen et 
al 0, 0, m IE 113 • The ballistic approximation allowed 
capturing essential behaviors such as front pulse accelera- 
tion, amplitude shrinking, and diminution of its duration. 
This lowest order model appeared accurate for small ta- 
pering factor, but may become insufficient when tapering 
becomes greater. Under this approximation, it was as- 
sumed that the front pulse adapts itself to a self similar 
solution at every impedance mismatch. This may clearly 
not always be true. Nevertheless, our data seemed to fol- 
low this scaling in a reasonable good manner indicating 
that the main parameters to characterize the quasisoli- 



7 



tary wave are the maximum force and the beads size, in 
addition to the mechanical properties of the beads. As a 
possible extension of this work, it would be interesting to 
extract a more complete behavior from a nonlinear wave 
equation that would take into account the tapering of 
the chain. In addition, comparisons of our experimental 
findings with numerical simulations are under way and 
should help to elucidate some of the hypothesis stated 



here. 
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